Large curvature perturbations near horizon crossing 
in single-field inflation models 
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We consider the examples of single-field inflation models predicting large amplitudes of the cur- 
vature perturbation power spectrum at relatively small scales. It is shown that in models with an 
inflationary potential of double-well type the peaks in the power spectrum, having, in maximum, 
the amplitude Viz ~ 0.1, can exist (if parameters of the potential are chosen appropriately). It is 
shown also that the spectrum amplitude of the same magnitude (at large k values) is predicted in 
the model with the running mass potential, if the positive running, n' , exists and is about 0.005 at 
cosmological scales. Estimates of the quantum diffusion effects during inflation in models with the 
running mass potential are given. 



PACS numbers: 98.80.Cq, 04.70.-s 



I. INTRODUCTION 



In last few months several papers appeared [l|, 0, H, 0] , 
in which single-field inflation models predicting (poten- 
tially) large amplitudes of the curvature perturbations 
on relatively small scales are discussed. It is shown in 
[l[ that large class of such models exists, namely, the 
models with a potential of hill-top type (the idea of the 
hill-top inflation was proposed, to author's knowledge, in 
the earlier work Q). In such models, the potential can 
be of concave-downward form at cosmological scales (in 
accordance with data) and be much flatter at the end of 
inflation when small scales leave horizon. Correspond- 
ingly, the amplitude of the perturbation power spectrum 
can be rather large. It is noticed in [l[ that the run- 
ning mass model, having the potential with the similar 
behavior, also can predict the large spectrum amplitude. 

Authors of [2 discuss also more general scenarios of 
producing large amplitudes of perturbation spectrum. 
They show the limitedness of the standard procedure of 
the potential reconstruction which can easily miss the 
potentials leading to large spectrum amplitude and to 
noticeable primordial black hole (PBH) production. 

In the recent paper [3[ it was shown that PBH produc- 
tion is possible in single-field models of two-stage type 
("chaotic + new"). The idea was proposed ten years ago 
in Authors of Q carried out the numerical calcula- 
tion of the power spectrum using the Coleman- Weinberg 
(CW) potential. 

In the present paper we continue a study of the prob- 
lems discussed in the previous works [l|, [H, Q . We inves- 
tigated thoroughly, as a particular example, the model 
of two-stage inflation with a potential of the double-well 
(DW) form, and showed that the characteristic features 
of the power spectrum in models of this type (such as an 
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amplitude and a position of the peak, a degree of tuning 
of parameters of the potential) are very sensitive to an 
exact form of the potential. Further, we carried out the 
numerical calculation of the power spectrum in a running 
mass model and showed that the spectrum amplitude at 
small scales can be rather large. Our calculation differs 
from the previous one Q in several aspects: we express 
the results through the values of parameters s, c, which 
are used nowadays and prove to be very convenient for a 
comparison with data; we studied, in details, the differ- 
ence in predictions of slow-roll and numerical approaches 
at high fc-values; we exactly specified the value of the 
positive running, n' , which corresponds to our spectrum 
prediction. In the final part of the work we investigated 
the quantum diffusion effects in a model with the running 
mass potential. 

A plan of the paper is as follows. In the second section 
we study predictions of two-stage inflation models with 
DW and CW potentials, with accent on a mechanism of 
the formation of peaks in the power spectrum. In the Sec. 
IIIII all aspects connected with an obtaining of the predic- 
tions of running mass inflation models are discussed. In 
the Sec. IIVI we present our main conclusions. 



II. 
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EXAMPLES OF THE POWER SPECTRUM 
WITH PEAKS 

A. Double-well potential 



This form of the inflaton potential having an unstable 
local maximum at the origin has been discussed many 
times in studies of eternal and new inflation. The main 
problem was to realize the initial condition for the new 
inflation when system starts from a top of the hill. Ten 
years ago the model of "chaotic new inflation" has been 
proposed Q , in which the system climbs on the top dur- 
ing dynamical evolution of the inflaton field with initial 
conditions coinciding with those of chaotic inflation mod- 
els. In the approach of Q the inflation has two stages, 
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FIG. 1: The solution of the background equation for inflation 
with the double- well potential ((TJ. The parameters of the 
potential are: v = 0.16286748m Pi , A = 1.7 x 1CT 13 . 
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FIG. 2: The time dependence of the parameter e and the 
combination 1 + e — 77 corresponding to the background field 
evolution shown in fig. [1] 



chaotic and new, and during transition from the first 
stage to the second the slow-roll conditions break down 
(in general). 

The potential has two parameters: 
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V(<t>) = T W-v J ) 



(1) 



The inflaton starts with the rather high value of <j> (we 



take (f>i n ~ 5mpi) and rolls down to the origin. The pa- 
rameter A is fixed by a normalization of the power spec- 
trum on experimental data, A ~ 10~ 13 . The evolution 
of the system strongly depends on the value of v: if v is 
finely tuned, <f> can spend some time near the origin, i.e. 
on the top, and then roll down to one of the two min- 
ima. In figs. [T]and[2k the time evolution for the inflaton 
and the parameter e for the definite values of the param- 
eters A , v are shown. One can see that, really, w at 
some period of time and, what is important, the slow-roll 
approximation is invalid (e ~ 1) just at the time of the 
transition from a rolling to a temporary stay at the top 
of the potential. 

It is well known that in situations when there is a fail- 
ure of the slow-roll evolution the perturbations on super- 
horizon scales can be amplified and specific features in 
the power spectrum can arise @, [i| EH EH, E2 (see also 
the recent paper [Hj]). In particular, in the earliest work 
where this problem was studied [8|, the inflation poten- 
tial with a sudden gradient discontinuity leading to the 
power spectrum of a step-like form was considered. All 
this means that the predictions of the slow-roll approxi- 
mation which are based on the assumption that pertur- 
bations reach an asymptotic regime outside the horizon 
cannot be trusted. 

The curvature perturbation on comoving hypersurfaces 
IZk, as a function of the conformal time r, is a solution 
of the differential equation (prime denotes the derivative 
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(<f) is the inflaton field). The standard initial condition 
for this equation, corresponding to the Bunch-Davies [l4| 
vacuum, is 



Ufc(r) = 



1 



—ikr 



2k 



aH < k 



(4) 



where u = zlZ. The variable u had been introduced in 

OElEa. 

The functions e and ij in eq. © are the Hubble slow- 
roll parameters defined by the expressions [TH 
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Outside the slow-roll limit these functions are not neces- 
sarily small. 



It had been demonstrated in that solutions of the 
equation ©, at k <^ aH, i.e., outside horizon, are well 
approximated by constant if the coefficient of the friction 
term, z' /z, doesn't change sign near the horizon crossing. 
In the opposite case, if z'/z changes sign at some time, 
the friction term becomes a negative driving term, and 
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one can expect strong effects on modes which left hori- 
zon near that time. In the present paper we study the 
corresponding features of the power spectrum, following 
closely the analysis of [ill ]. 

According to eq. ([3]), z' / z is proportional to 1 + e — r\ 
and the comoving Hubble wave number aH. The time 
dependences of these functions are shown in fig. (2b- One 
can see that the interruption of inflation correlates with 
the change of the sign of 1 + e — r\. 

The time evolution of curvature perturbations for sev- 
eral modes is shown in fig. [31 It is clearly seen that the 
perturbations TZk for different modes freeze out at dif- 
ferent amplitudes. The mode which crosses horizon near 
the moment of time when the sign of 1 + e — rj changes 
(i.e., near t ss l.bmpi) freezes at maximum amplitude, 
due to the exponentially growing driving term in eq. 
(which is most effective just for this mode). It leads to 
the characteristic peak in the power spectrum "PTj(fc), 

Vn{k) = j^\n k \ 2 , (6) 

shown in fig. 21 

The calculations of TZk (fig- G2) are carried out up to 
the end of inflation, and the power spectrum in fig. [4] 
also corresponds to this moment of time. We estimate 
approximately the reheat temperature in our case as ~ 
(Aw 4 ) 1 / 4 - 10 14 GeV. The horizon mass at the be gmnmg 
of radiation era is 

M M ~ 10 "g(l^X) 2 ~10» g , (7) 

and maximum wave number, which equals the Hubble 
radius at the end of inflation, is 

(M \ 1/2 

fcead = Geq-ffeq ( jf 1 ) ~ 10 23 Mpc" 1 . (8) 

In the fig. [S]we show the power spectrum calculated for 
two stages of its evolution: at horizon exit (HE) and at 
the end of inflation (END). In the same figure the result 
of the calculation with an use of the slow-roll formulae is 
also shown. The peak of HE curve at the region of large k 
is due to a failure (for TZk) to reach the asymptotic limit. 
One can see also that the slow-roll approximation is too 
crude to describe perturbations at the end of inflation 
(this our conclusion agrees with general statements of 

m 

B. Coleman- Weinberg potential 

The Coleman- Weinberg potential has the form [l9j : 

It looks very similar to the previous one, but the impor- 
tant difference is its behavior near the origin. Namely, 
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FIG. 3: A time evolution of the curvature perturbation TZk if) 
for several different values of wave number k during inflation 
with the DW potential. The parameters of the potential are 
the same as in fig. [l] 
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FIG. 4: The numerically calculated power spectrum Vn(k) for 
the model with the potential ([1]). Parameters of the potential 
used in the calculation are the same as in fig. [T] 
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FIG. 5: The power spectrum Vnik) for the model with po- 
tential (jTJ) . END: V-nik) is calculated at the end of inflation; 
HE: Vn{k) is calculated at the time of the horizon exit; SR: 
the slow-roll result. 
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FIG. 6: The result for the power spectrum Vn(k) calculation 
for the CW potential ([9|, for two sets of parameters. Left 
peak is for v = 1.113Mp, A = 5.5 x 10~ 13 . For the right peak, 
v = 1.112Mf, A= 2.4 x 10" 13 . 



the CW potential behaves as A + _B</> 4 ln(0/i>) near the 
origin, i.e., it is more flat near zero, in comparison with 
the DW potential. Therefore, it has more e- folds of "new 
inflation" 6] and, as a consequence, the peaks of the 
power spectrum (arising, as in the previous case, due to 
the temporary interruption of inflation) correspond to 
relatively smaller k values. Besides, at the beginning of 
new inflation when the very flat region of the potential 
near zero is crossed by </>, quantum fluctuations with a 
particle creation can be large (e.g., see below, Sec. IIII C[) 
and must be taken into account. 

In fig. [6] two examples of the power spectrum calcula- 
tions are shown for two different sets of parameter values. 
As before, the peaks are very distinct, although their am- 
plitudes are smaller. 



C. Possibilities of PBH production 

One can see, in particular, from fig. [H that, in princi- 
ple, the production of primordial black holes (about these 
objects, see original works [2(1 HH and reviews [22l |23|) 
can be rather large in single-field inflation models with 
potentials of double-well type. The main conclusion is 
that it requires rather large fine tuning of parameters of 
the potential. The characteristic PBH mass is estimated 
by 
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BH 



M h = M, 



hi 



^cnd 
^pcak 



(10) 



and in the case of the spectrum of fig. [4] A/bh ~ 10 7 g. 
In the CW case, fig. M BH ~ 1OOM for the left peak 
and Mbh ~ 10 27 g for the right peak (but the amplitudes 
of the spectra are too small). 

Recently, it has been shown Q that inflation with 
CW potential is capable to produce significant number of 
PBHs: the parameter v can be chosen (by finest tuning!) 



in such a way that inflaton field makes several oscillations 
from one minimum to another before it climbs on the top 
and "new inflation" starts. 

In the present paper we do not consider possibilities of 
a constraining of the peak amplitudes, based, for exam- 
ple, on the effects of PBH evaporation in early Universe 
(see, e.g., 0). 



III. 



RUNNING MASS MODEL 



A. Main assumptions and approximations 

We consider in more detail a case of the running mass 
inflation model [H, H, [27|, [H, 1<1 [H HI which predicts a 
spectral index with rather strong scale dependence. The 
potential in this case takes into account quantum correc- 
tions in the context of softly broken global supersymme- 
try and is given by the formula 



V 



Vo + ^m 2 (ln 



(11) 



The dependence of the inflaton mass on the renormaliza- 
tion scale cf> is determined by the solution of the renor- 
malization group equation (RGE). 

1. The inflationary potential in supergravity theory is 
of the order of M^ f , where Mj n f is the scale of super- 
symmetry breaking during inflation. In turn, the mass- 
squared of the inflaton (and any other scalar field) in 
supergravity has, in general, the order of the square of 
Hubble expansion rate during inflation, 



Vo 
3M£' 



(12) 



We suppose, for simplicity (see [251, |2 
Mi n f ~ M s , where M s is 4-1 1 



breaking in the vacuum, 



_ .. _„ M Si), that 
the scale of supersymmetry 



M s 



y/fhsMp ~ 10 n GeV - 3 x KT 8 M P (13) 



(rh s is the scale of squarks and slepton masses, m s ~ 3 
TeV). These assumptions give the scale of the inflationary 
potential: 



Vo - M s 4 - 10 _30 M 4 



H, 



1Q- 15 M P . 



(14) 



2. RGE for the inflaton mass is the following (we 
consider a model [H, [29[ of hybrid inflation using the 
softly broken SUSY with gauge group SU(N) and small 
Yukawa coupling): 



m 2 (t) = to 2 , — Attiq 
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TOq and rfiQ are, correspondingly, the inflaton and gaugino 
masses at <j> = Mp, 
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Ba 
2tt ' 



(16) 
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a is the gauge coupling constant, a — .g 2 /47r. A and 
B are positive numbers of order 1, which are different 
for different variants of the model, even if they are based 
on the same supersymmetric gauge group SU(N) (it de- 
pends on a form of the superpotential, particle content 
of supermultiplets, etc). We use in the present parcr the 
variant of [29] and, correspondingly, put everywhere be- 
low A = 2 and B = N = 2. 

3. A truncated Taylor expansion of the potential 
around the particular scale tf>o (in our case, O is the 
inflaton value at the epoch of horizon exit for the pivot 
scale k 0.002/i Mpc" 1 ) is 

(ln(0 o )) - c-rrjL In — 
Mp 0o 

Here, constant c is defined by the equation 
Vb dm 2 



v{4>) = v + ^ 



(17) 



M 2 



(18) 



0=00 



In turn, a Taylor expansion of eq. (|15j) up to linear 
terms gives (to = In -^f-)' 



m 2 (t) = m 2 (t ) 
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From eqs. (]T8[) and (JT9J) we obtain the expression for the 
constant c, 
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FIG. 7: a) Evolution of the inflaton field <p(hxa) in the run- 
ning mass model, b) The dependence of the parameter 77 
on a value of the field (f>. For both plots, Hj = 10 -15 Mp, 
c = 0.062, s = 0.040. 
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(21) 



It appears (see fig. 03) that in our example <po ~ 
IQ- 10 M P , so, t - lnlO" 10 - (-23). Assuming that 
a ~ 1/24 (as in SUSY-GUT models), one has a ~ 
^7^4. In such a case, c ~ 4<5 ~ 0.06. 

If we would keep terms of higher order in t — to = In 4^ 
in the Taylor expansion of m 2 (t) in eq. (fH)|) we would 
see that the real expansion parameter is do In rather 
than In . The smallest value of </>, </> cn d 

~ 10 _ 16 Mp (see fig. [7b). Even for such value of 
the expansion parameter is rather small, 



4>o 

in our case is 

^end; 



do In 
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So In 10 



(-0.1) 



(22) 



Having this in mind, we will use the linear approximation 
for the inflaton mass (eq. §19\i ) in the entire region of 
inflaton field values exploited in the present paper. 

Following the previous papers, we introduce also an- 
other parameter, 



where is the inflaton value corresponding to a maxi- 
mum of the potential. This parameter connects the field 
value (f>o with the Hubble parameter during inflation and 
with the normalization of the CMB power spectrum: 
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4. The minimum value of the inflaton field which cor- 
responds to the end of inflation can be determined from 
the approximate equation [2i| 
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Using RGE, one obtains from this formula the relation 
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Substituting here A — 2, ml = |toq| = Vo/iMp, one has 
finally the approximate expression for </> e nd) 



s = cln I — 



(23) 



"end 



Mp exp 



ao V \/3 



(27) 
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which shows that the minimum field value is very sensi- 
tive to the value of the model parameter a and, in our 
case, does not depend on Vq- More exactly, the condition 
X] = 1 means the end of the slow-roll part of inflation. We 
suppose, as usual (see, e.g. [25|, (26|) that in reality in- 
flation ends by hybrid mechanism, and the critical value 
of inflaton field, cr , is determined by the value of the 
Yukawa coupling A (in spite of the inequality A 2 <C a). 
One can check [29| that the value of A can always be cho- 
sen such that 4> CT < (j> cn d and slow-roll ends before the 
reaching of <p CI . 

One should note that the accuracy of the approximate 
formula (f27|) is not very good. Luckily, in the approach 
based on the numerical calculation of the power spec- 
trum there is no need to use it, because the value of </> C nd 
appears in a course of the calculation (fig. [7b) . 



B. Power spectrum of curvature perturbations 



t data [32j, |33j , includ- 
, leads to the following 



An analysis of CMB anisotro 
ing other types of observation [3 
main qualitative conclusions: 

i) the power spectrum of scalar curvature perturba- 
tions is red, i.e., the spectral index is negative, 



n = 0.97 ±0.01 



(28) 



ii) observations are consistent, or, at least, are not in 
contradiction with the small positive running of the spec- 
tral index, n < 0.01; 

Hi) the contribution of tensor perturbations in the 
value of the spectral index is small (< 10 -2 ) and, as a 
result, n ~ 1 + 2rj ; it means that r\ is negative, and 
the potential must be concave-downward (i.e., of hill-top 
type), while cosmological scales cross horizon during in- 
flation 

One should note that, strictly speaking, the conclusion 
Hi) is not grounded firmly enough. According to the 
recent analysis [35] , the present data still admit any sign 
of r\ and V" . 

These conclusions constrain the possible values of the 
parameters s and c. Approximately, for cosmological 
scale one has 



no — 1 2(s — c) 



2sc . 



(29) 



From the conclusion Hi) it follows that c > (it is con- 
sistent with eq. (JM])), from the positivity of n' (the con- 
clusion ii) ) it follows that s > 0. At last, the conclusion 
i) leads to the inequality s < c. 

We choose for the power spectrum calculation the fol- 
lowing values: 



0.062 , s = 0.040. 



(30) 



These numbers correspond, at cosmological scales, to the 
following values of slow-roll parameters: 
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c~ (-0.02) 



(31) 



a, 

so 



lg k [MpcT 1 ] 

FIG. 8: Power spectrum Vn{k) in the running mass model, 
calculated numerically (solid line), by the approximate ana- 
lytic formula (l32jl (long-dashed line) and using the Stewart- 
Lyth extended slow-roll approximation (short-dashed line). 
The parameters of the potential are the same as used in fig. 
[7J The arrow shows the value of fc cn d- 



that seems to be consistent with the present data 

To check the validity of the slow-roll approximation, 
we calculate the spectrum by the three ways: i) using 
the approximate analytic slow-roll formula 



= exp 



2s / 



zAN(k) 



1) - 2cAN(k) 



(32) 



(AN(k) = ln(fc/fco); this expression is easily derived from 
the simplest slow-roll prediction 



H 2 



'pi 



(33) 



aH=k 



which gives the power s pec trum to leading order in the 
slow- roll approximation [361]): ii) using the Stewart-Lyth 
approximation [371 ] . which is valid to first order in the 
slow-roll approximation, 



V](\k) = [1 



27T U 



(34) 



iH=k 



c 



-0.73 



Hi) by numerical integration of the differential equation 
for K k , eq. ©• 

The results of the calculations are presented in figs. El- 
El Fig. [7] shows the evolution of the inflaton field <j> with 
the scale factor and a growth of the slow-roll parameter r\ 
with a decrease of 4> from (j)Q to 4> e nd- The power spectrum 
is shown in fig. [8] for a broad interval of comoving wave 
numbers. It is clearly seen that near the end of inflation, 
when 

4> - 10~ 16 M P , k - fc cnd = flendffend = SxlO^Mpc" 1 , 

the slow-roll formulae are inaccurate: they strongly un- 
derestimate values of V-r. To illustrate this point more 
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FIG. 9: Dependence of Vn{k,a) calculated numerically and 
Viz from the Stewart-Lyth formula. The comoving wave num- 
ber for this figure is k = 10 15,8 Mpc -1 . 



the coefficient H 2 ^ can be rewritten as y 4 ^/ At , and the 

evolution of tp on timescales At > H~ x can be described 
by a finite-difference form of the eq. (|36j) , 

<p(t+At)-<p(t) = -^V\y)M+^yfWKtm- (38) 

The condition for the deterministic evolution is (see, e.g., 
0) 



— |V'foO|At»— Vff^At , At = H~\ (39) 

Using the slow-roll connection between V and 7?, one 
obtains 

i^)i>>|^ 3 = ^ 3/ >)- (40) 



clearly, at the next figure we show the comparison of two 
curves: aii-dependence of the spectrum calculated nu- 
merically for a definite value of k, k = 10 15 8 Mpc -1 , and 
aH- dependence of the Stewart-Lyth spectrum. It is seen 
that the numerical spectrum at the moment of crossing 
horizon (when aH — k) is already almost asymptotical, 
and its value distinctly exceeds the corresponding value 
predicted by Stewart-Lyth formula. 



In the approach of [43j the coarse-grained field is con- 
sidered as a perturbation of the classical solution tp c \ 
(which is the solution of the Langevin equation without 
the noise), 



tp(t) = tp cl (t) + S(p!(t) + 6cp 2 (t) + .. 



(41) 



Here, the term Sipi(t) depends on the noise at the power i. 
It is assumed that the Hubble parameter in the Langevin 
equation depends only on the coarse-grained field tp, 



C. Quantum diffusion effects 



H 2 (tp) = 



1 



3MJ, 



V(tp). 



(42) 



We calculated the curvature perturbations in terms of 
the classical trajectories of a scalar field associating, in 
particular, points in a field space with definite numbers 
of e-folds from the end of inflation. This description 
becomes incorrect if the quantum diffusion destroys the 
classical evolution of the field. In this case we should use 
the methods of stochastic inflation. The latter approach 
operates with the coarse-grained field, which is defined to 
be spatial average of the field 4> over a physical volume 
with size larger than the Hubble radius H . 

In the slow-roll approximation, the evolution of the 
coarse-grained field tp is gov erned by the first order 
Langevin-like equation [38l l39l. lifjl. IlO] 



Correspondingly, the Hubble parameter can be expanded 
perturbatively, 



H{tp) = H cl + H^Stpi + 5<p 2 ) + ^-Stpf 



H(tp cl ) 



3Mf, 



(43) 



(44) 



This approach permits to calculate the mean value of 
the total number of e-folds, (N) , and to compare it with 
the corresponding "classical" number, 



1 H 3 / 2 
V+^-V'(tp) = ^—^(t), 



3H 



2tt 
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Vend 



NSi = -- 



1 f , H cl 



(37) 



Here, £(t) is a random noise field, and angular brackets 
mean ensemble average. The term -^V'(tp) describes 
the deterministic evolution of the field tp, in the absence 
of the noise term H 2 ' £(t). The solution of eq. (|3"6"|) , in 
the absence of the noise term, is the deterministic slow- 
roll trajectory tp SI (t). Going to finite time differences, 
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dtp cl . (46) 
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FIG. 10: a) Comparison of V and V 3/2 for the running 
mass model (see eq. ((40])). b) The same figure in different 
scale: variable x is connected to the field value by the relation 
^=(1-10-)^. 
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FIG. 11: The result for the calculation of 8tp/tp for different 
values of ipi n — <p*(l — 10~ x ). From bottom to top, x equals 
3, 4, 5. Here, Sip = y/ (5(p\) + (Sip 2 ). 



Besides, one can calculate the mean value of the Gaussian 
probability distribution function for the coarse-grained 
field and to see how it behaves as a function of the current 
field value. 

The interval of inflaton field values for which the in- 
equality holds and, therefore, the deterministic evo- 
lution dominates, is shown in fig. UOa .b. The variable x 



in fig. [10b is defined by the relation 

^ = 1 - 10~ x (47) 

(<y5* is, as before, the point of a maximum of the potential 
V(ip)). It is seen from the figure that the constraints on 
the inflaton field following from the condition (T4"0|) are 
not too severe, 

3 x 10- 19 M P < ip < tp* - 10"V ■ (48) 

An accuracy of the perturbative expansion (|41|) for 
the monomial potential of stochastic inflation models was 
studied in [4J| . Here we study this accuracy for the run- 
ning mass potential, using the parameters Vq, s, c, (f C nd 
introduced above. The results of the calculation of (Sipi) 
and {5ip2) are shown in fig. [TlJ As one can see, the per- 
turbative expansion is good if the starting value of the in- 
flaton field, is chosen to be not too close to the value 
of <p at the maximum of the potential. More exactly, the 
parameter x, defined in eq. (|47p . must be smaller than 
4-^4.5. The starting value <p- lu corresponds to a beginning 
of the evolution, i.e., <5y>i(ii n ) = ^^(^in) = 0. 

If the evolution is deterministic, the correction to a 
classical e-fold number A^ 1 , given by eq. (|46|) , is small. 
To estimate analytically the upper limit for this correc- 
tion, we used the analytic expressions for (8<p\) and {Sip?) 
derived in [43| . keeping in them leading terms only. Ac- 
cording to these expressions, the following inequalities 
hold: 

( ^ 2> < (^) (^) m%^ ~ 

< lQ- w+2x <p , (49) 

(Spl) < (^V^- ~ 10-3°+^. (50) 

Using these upper limits, one can estimate the corre- 
sponding upper limits of two integrals in the expression 
for SN T (|3B|), The result is the following: 

^ J {S^Wox < lQ- 30+2x ; (51) 

fend 

fend 

It is clear from the inequalities (|52|) that the quantum 
correction to e-fold number, SNt, is quantitatively small 
even if the value of x is as large as 10. But only if x < 
4 -;- 4.5 , and the perturbative expansion, eq. (I4ip . is 
valid, one really can be sure that 

8N T < , (53) 
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FIG. 12: Calculation of stochastic effects for various val- 
ues of initial field ifi n . a) Thin dashed curve: <fi n = (1 — 
1.5 x lO -5 )^*, no volume effects included; thick dashed curve: 
(fin = (1 — 1.5 x 10~ 5 )y*, volume effects included; solid thick 
curve: tp- m = (1 — 3 x 10~ 5 )</p* , volume effects included, b) 
Thick long-dashed and thick short-dashed curves correspond 
to the cases with and without inclusion of volume effects, re- 
spectively, for if i n = (l — 3 x 10~ J )</p». 

and the evolution is deterministic. 

However, this analysis is still not complete: one must 
check also the position of the mean value of the probabil- 
ity distribution function for the coarse-grained field ■ 
The calculation, with taking into account the volume ef- 
fects, leads to the results shown in fig. [l"2"kb. It is seen 
from the figure that in this case also, as in the calcula- 
tion of the e-fold number correction, the correct choice 
of the initial condition plays a decisive role: there are no 
effects of a "walk" of the mean value around the maxi- 
mum of the potential (such effects were noticed in (43[) if 
evolution starts from the point which is far enough from 
the maximum (x < 4.5). Of course, the realization of 
the initial condition of this kind is, in itself, a problem. 
Supposedly, it could be provided by the previous history 
of eternal inflation [25|, [26| . 

IV. CONCLUSIONS 

1. It is shown, by numerical methods, that in the 



single-field inflationary model with a simple double-well 
potential the parameter values of this potential can be 
chosen in such a way that the power spectrum of curva- 
ture perturbations Vn(k) has a huge peak (with ampli- 
tude ~ 0.1) at large k (and the right normalization and 
monotonic behavior at cosmological scales). The peak 
arises due to temporary interruption of the slow-roll near 
points of the minimum of the potential, ±u. The corre- 
sponding mass of PBHs produced in early universe in a 
case of the realization of such power spectrum is about 
10 7 g. 

The analogous behavior of the power spectrum was ob- 
tained by authors of [3] in a model with CW potential. 
There are some important differences in the results of 
Q and ours, in the peak amplitude and PBH mass, con- 
nected, in particular, with a large flatness near the origin 
in a case of the CW potential. 

2. It is shown that the inflation model with running 
mass potential predicts a rather large amplitude of the 
power spectrum of curvature perturbations (~ 0.1) at k- 
values ~ 10 16 Mpc -1 . For such a prediction, a very small 
positive spectral index running at cosmological scales is 
necessary, n' ~ 0.005, as well as a small negative value 
for the slow-roll parameter r\ (ps —0.02). Both this num- 
bers do not contradict with data. It is shown also that 
for an obtaining the correct quantitative results for the 
power spectrum at largest fc-values an use of numerical 
methods is required because, in general, slow-roll formu- 
las are not accurate enough at the end of inflation, where 
r\ « 1. 

3. Quantum diffusion effects in a model with the run- 
ning mass potential are studied in details. It is shown 
that inflationary evolution of the universe in a model 
with a scalar field and the running mass potential can 
be described by the classic deterministic equations, and 
for a possibility of such a description the correct choice 
of the initial conditions is crucial. Concretely, an initial 
value of the inflaton field (at the beginning of the evolu- 
tion) should not be too close to a point of the maximum 
of the potential. If this condition is satisfied, the quan- 
tum corrections to a total e-fold number and to a position 
of the mean value of the probability distribution function 
are small. 
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